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Abstract 

In this report we investigate the multiscale hydrodynamical response of a liquid as a function of mixture 
composition. This is done via a series of molecular dynamics simulations where the wave vector dependent 
viscosity kernel is computed for three mixtures each with 7-15 different compositions. We observe that the 
nonlocal viscosity kernel is dependent on composition for simple atomic mixtures for all the wave vectors 
studied here, however, for a model polymer melt mixture the kernel is independent of composition for large 
wave vectors. The deviation from ideal mixing is also studied. Here it is shown that a Lennard- Jones 
mixture follows the ideal mixing rule surprisingly well for a large range of wave vectors, whereas for both 
the Kob- Andersen mixture and the polymer melt large deviations are found. Furthermore, for the polymer 
melt the deviation is wave vector dependent such that there exists a critical length scale at which the ideal 
mixing goes from under-estimating to over-estimating the viscosity. 
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devices. Recent studies 



Kob-Andersen Lennard- Jones 

caa/^aa 1 1 

(-bb/€aa 1/2 1/2 

£ab/ £aa 0.8 yi72 

o'aa/o'AA 1 1 

obbIoaa 0.88 0.88 

& abI v aa 0-8 0-94 

TABLE I: List of the Kob-Andersen and Lennard- Jones parameters used in this work. 

Hydrodynamics on very small length scales has become an important research area because 
it is believed to hold the Jcey^ to understand the many new phenomena observed in nanofluidic 

have shown that the spatial correlations in the fluid reduce the 
shear stress compared to the stress predicted via a local response function. The nonlocal response 
is described via generalized hydrodynamics [3], in which the response function is a wave vector 
dependent quantity. The wave vector dependent viscosity, i.e. the viscosity kernel, accounts for 
the momentum flux due to nonzero strain rate. It has been found for one component fluids through 
molecular dynamics simulations and it was found that it follows a simple functional form reasonably 
well for atomic, diatomic and polymer fluids [3H6|]. By now the effect of the multiscale response is 
understood fairly well for a range of simple single component fluids 

In order to investigate the multiscale hydrodynamical response as a function of composition, 
we here evaluate the viscosity kernel for three different two component mixtures using molecular 
dynamics. In details, these are (i) a Kob-Andersen (KA) mixture {{J, (ii) a Lennard- Jones (LJ ) 
mixture using the Lorentz-Berthelot mixing rule, and (hi) a model polymer melt mixture 1 1 Of ] . 
In all the simulations the particles interact through the Lennard-Jones cut and shifted potential 
Uuifij) = 4e (cr/rij) 12 — (<r/rjj) 6 — U(r c ) for rij < r c , where is the distance between particle 
i and j, a is a length scale, e is an energy scale, r c is the interaction range (cut-off) and U(r c ) 
is the unshifted potential at r c . The values of a and e are different for the KA and LJ mixtures 
depending on the pair of particles that interact, see Table HI 



t of the multiscale respc 
[3-6] and glasses 



melt. The latter is also referred to as the Weeks-Chandler- Andersen pair potential [11]. In the 
polymer melt the particles (or beads) are bonded via the Finite Extensible Nonlinear Elastic 
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(FENE) potential [l0(] U FE ne = -kR ln 1 - (nj/Ro) 2 /2, where k = 30e/a 2 and R = 1.5a. 
The polymer melt is composed of two types of molecules, one with two beads, component B, and 
one with ten beads, component A. In what follows we give all quantities in terms of Lennard- 
Jones reduced units, for example, reduced distance r*- = r^/cr and number density p* = pa 3 . For 
simplicity of notation, we will hereafter omit the asterisk. 

The simulations are carried out at an average reduced pressure of p = 1 and temperature 
T = 2.5. The target pressure was obtained via an anisotropic Berendsen barostat [12] such the 
simulation box was extended in the x direction only. The temperature was controlled via a Nose- 
Hoover thermostat [ijl, []J ] • 

The expression for the wave vector dependent viscosity can be found from the generalized Navier- 
Stokes equation and is given in terms of the transverse momentum current density autocorrelation 



function C± (k, t) [la. Il6l] 



. C ± (k,t = 0)-iuC ± (k,u}) 

mk,uj) = , (1) 

V ' C ± (k,co)k 2 /p V ' 

where p is the mass density, k is the z component the wave vector, i.e. 

k = 2-Kn/L z , n= 1,2, ... , (2) 

where L z is the simulation box length in the z direction, and C±(k,t) = { J y (k, 0) J y (k, t))/V. The 
transverse momentum density is here defined via J y (k,t) = miV y (t)e tkz , where rrii and v y are 
the center of mass and center of mass velocity of molecule or particle i. We note that C±(k,u)) 
is the Fourier-Laplace transform of C±(k,t), that is, C±(k,ui) = f£° C±(k,t)e lult dt. Also note, 
because the barostat is anisotropic and only varies the simulation box in the x direction L z is 
constant in Eq.([2]). 

In Fig. Q] we have plotted the viscosity kernel data in the limit oj — > for different composition 
fractions of A, xa = NA_/Nt where iV^ is the number of A particles and Nt is the total number of 
particles. 



Recall, that in the polymer melt we denote the 10 bead polymer A. From Fig. [I] it is observed that 
for the KA and LJ mixtures the hydrodynamical response is dependent on the exact composition 
for all length scales studied here. This is not the case for the polymer mixture: here the fluid 
response is largely independent of the fluid composition for sufficiently large wave vectors. In order 
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°i i * 6 rtr 




°A 1 i 3 — -4 6 6 1 



FIG. 1: Viscosity kernel data (symbols) for the Kob- Andersen mixture (a), the Lennard- Jones mixture (b) 
and the polymer mixture (c). In (a) and (b): xa — (filled squares), 0.8 (upwards pointing triangles), and 
1 (downwards pointing triangles). In (c): xa = (filled squares), 0.25 (upwards pointing triangles) and 1 
(downwards pointing triangles). The error bars are the standard error. Lines represent the best fit of the 
data to Eq. ([3]). The inset windows depict the kernels for large k. 
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to decrease the statistical error in the further analysis we fit the data to a Lorentzian functional 
form 

where a and f3 are fitting parameters that depend on the fraction of A, xa, and t]q(xa) is the zero 



frequency viscosity for k = 0. Eq. ([3]) has been s 
of different non-charged single component fluids 



lown to fit viscosity kernel data well for a range 
(J. The result of the fitting is also depicted in 



Fig. Q] where the inset windows show the data and the fits for large values of A;. It is seen that Eq. 
([3]) fits data well for all three systems, however, we stress that the agreement is not satisfactory 
for large values of k in the cases of Kob- Andersen and Lennard-Jones mixtures. With this in mind 
we, will from now on use the fitted values of the viscosity kernels rather than the raw MD data. 
We can define the k dependent excess viscosity as 

rj E (k, x A ) = rj(k, x A ) - r] ld (k, x A ) , (4) 

where rj l(i (k, xa) is the ideal part of the viscosity kernel given by an Arrhenius type mixing rule 13, 

Q 

7 ] id (k,x A ) = r ] A(k) XA ilB(k) 1 - XA (5) 

Here r)A{k) and rjsik) are the viscosity kernels of pure A and B, respectively. The excess viscosity 
has been fitted to various simple mixing models including the Kendall-Monroe model |l8i], the 



^ederer model, see for example Ref. [19|], and an extended version of the Grundberg-Nissan model 
2o| . We have found that the fourth order McAllister model 21 1 fitted the data best. This model 



is originally written in terms of the dynamical viscosity v = r]/ ' p and may readily be extended to 
include the wave vector dependency, that is, 

\nv{k,x A ) = x\ln[iy A (k)] + A x 3 AX B ln[M 31 (k)}+6x 2 AX 2 B ln[M 2 2(k)} + 

\xax\ ln[Afi 3 (fe)] + x% ln[i/ B (k)) - \u(xa + x B m r ) + 
„3„ „w, /3 + m r \ 2 2 fl + m r 



Ax A x B In [-^— ) + ^A^B ^ 1^^— ) + 
4x A x B In ( 1 + ^ m " r + X A B ln(m r ) , (6) 

where m r is given by tub / tha and M31 , M22 and M13 are the wave vector dependent McAllister 
coefficients. Recall, since we study binary mixtures xb = 1 — xa- From Eq. one can easily 
extract the McAllister excess kinematic viscosity using M31 , M22 and M13 as fitting parameters 
Rather than comparing the absolute excess viscosities, we compare the relative deviation from ideal 
mixing using rj E (k, XA)/r](k, xa) = 1 — rj ld (k, XA)/r]{k, xa)- This is done in Fig. [2j 



We also note that the normalized excess viscosity is a measure of the relative deviation from ideal 
mixing. Firstly, it is observed that the KA and polymer mixtures feature large deviations from the 
ideal mixing rule compared to the LJ mixture, that is, the Lorentz-Berthelot mixing rule features 
a more Arrhenius-like behavior compared to the KA mixing rule for all wave lengths. Secondly, for 
the KA and LJ mixtures the relative excess viscosity varies very little with respect to wave vector 
meaning that the relative difference between the ideal mixing response function and the actual 
response is weakly wave vector dependent. For the polymer mixture, Fig. c), we observe that 
the relative excess viscosity is strongly wave vector dependent. This is also indicated in Fig. Q] c) 
where it was observed that the kernel is independent of composition for large k, but not for small 
k. This complex wave vector dependency means that there exists a critical length scale l c at which 
the ideal mixing term goes from under-estimating to over-estimating the viscosity. For the melt 
studied here it happens around l c = 1.7. This behavior indicates that for the polymer melt the 
local viscous response is independent of chain length at short wave lengths. 

To discuss this special behavior further we have plotted the McAllister coefficients as functions 
of wave vector in Fig. [3l 

We see that Mij(k)/Mij{0), {ij} = {31}, {22} and {13}, fall on a master curve in the case of KA 
and LJ mixtures. This means that the three functions, Mij(k), may be described by a single master 
function which is directly proportional to any of the three Mij(k). For the polymer mixture, on 
the other hand, this is not the case. Note, that all the McAllister coefficients follow a Lorentzian 
form, see Eq. (|3|), but for the polymer mixture the parameters a and /3 are dependent on the index 

In this report we have investigated the multiscale hydro dynamical viscous response as a function 
of fluid composition. This was done through the viscosity kernel that was computed via equilibrium 
molecular dynamics simulations. We studied three different mixtures, namely, (i)a Kob-Andersen 
mixture, (ii) a Lennard-Jones mixture using Lorentz-Berthelot mixing rule, and (hi) a polymer 
melt mixture. We observed that the viscosity kernel is independent of the wave vector for large 
wave vectors in the case of a polymer melt, that is, the hydrodynamical response at these length 
scales is independent of the composition. This was not the case for the simple Kob-Andersen and 
Lennard-Jones mixtures. The deviation from ideal mixing is low in the case of the Lennard-Jones 
mixture, i.e. the Lorentz-Berthelot mixing rule agrees reasonably well with the predictions from 
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ideal mixing for all wave vectors studied here. This is not the case for Kob- Andersen and polymer 
mixtures. Finally, the relative deviation from ideal mixing is relatively wave vector independent 
in the case of the Kob- Andersen and Lennard-Jones mixtures. However, for the polymer mixture 
this deviation shows a strong wave vector dependency, since the ideal mixing rule does not predict 
the largely wave vector independent behavior of the kernel, Fig 1 c). 
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FIG. 2: Relative excess viscosity as function of composition and wave number, n. a) The Kob- Andersen 
mixture, b) the Lennard- Jones mixture and c) the polymer mixture. Full li nes are the corresponding fits 
using Eqs. (O and (JSJ. 
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FIG. 3: McAllister coefficients as function of wave vector, k. Filled squares represent the results for 
the Kob- Andersen mixture, upward pointing triangles the Lennard- Jones mixture and downward pointing 
triangles the polymer mixture. Lines serve as a guide to the eye. 
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